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We describe cosmic gamma-ray and neutrino signals of dark matter annihilation, explaining how 
the complementarity of these signals provides additional information that, if observable, can en- 
lighten the particle nature of dark matter. This is discussed in the context of exploiting the separate 
galactic and extragalactic components of the signal, using the spherical halo model distribution of 
dark matter. We motivate the discussion with supersymmetric extensions of the standard model 
of particle physics. We consider the minimal supersymmetric standard model (MSSM) where both 
neutrinos and gamma-rays are produced from annihilations. We also consider a gauged B — L, 
baryon number minus lepton number, extension of the MSSM, where annihilation can be purely 
to heavy right-handed neutrinos. We compare the galactic and extragalactic components of these 
signals, and conclude that it is not yet clear which may dominate when looking out of the galactic 
plane. To answer this question, we must have an understanding of the contribution of halo substruc- 
ture to the annihilation signals. We find that different theories with indistinguishable gamma-ray 
signals can be distinguished in the neutrino signal. Gamma-ray annihilation signals are difficult to 
observe from the galactic center, due to abundant astrophysical sources; but annihilation neutrinos 
from there would not be so hidden, if they can be observed over the atmospheric neutrinos produced 
by cosmic rays. 
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I. INTRODUCTION 

Understanding the particle nature of dark matter is a major problem of modern physics. While indirect gravitational 
evidence of its presence is plentiful, unambiguous identification of its particle properties is actively sought. Presumably, 
this population of particles is a relic that was produced spontaneously in, and is left over from, the big bang. It is 
common for extensions of the standard model that contain a dark matter candidate to have a Z2 symmetry that 
stabilizes the particle from decay, but does not prevent the particle from self-annihilating. Often in viable theories, 
the annihilations in the early Universe are important to bring the dark matter's relic density down to the observed 
density today. Once the rate of expansion of the Universe becomes larger than the dark matter annihilation rate, 
annihilations become rare, and the dark matter relic is said to freeze out at this time. However, if dark matter 
annihilates, then rare annihilations continue to occur today, predominantly in the densest regions of the Universe. 
Unambiguous identification of cosmic radiation from these annihilations would not only provide valuable information 
about the particle nature of dark matter, but also information about the distribution of the matter responsible for 
the signal. 

The information about dark matter from indirect detection of its cosmic annihilation radiation is complementary to 
that gleaned from other current experiments. Indirect detection experiments constrain the dark matter particle mass, 
annihilation cross section, and annihilation spectrum. Meanwhile, direct detection experiments attempt to observe 
dark matter-nucleon interactions in a laboratory detector, and constrain the dark matter mass and its nucleon- 
scattering cross section. Particle accelerators try to detect dark matter production from particle collisions, where the 
dark matter would be manifested in missing transverse energy. Again, events of this kind provide information on the 
dark matter's particle mass, but also on the processes that led to the creation of the dark matter. The results from 
these experiments provide different constraints for particle physics models, and provide consistency checks for one 
another. 

Today, indirect detection experiments are looking for high-energy cosmic rays, gamma rays, and neutrinos produced 
from dark matter annihilation. The propagation of produced cosmic rays within the galaxy is difficult to describe pre- 
cisely, complicating the prediction of the signal's properties. Also, energy losses do not allow extragalactic cosmic rays 
to reach us. However, neutrinos and gamma-rays can be observed from their sources, both galactic and extragalactic. 
This allows their observed signals to probe not only the particle physics that produced them, but also the distribution 
of their sources: the matter in our galaxy or extragalactic large scale structure. A sample of current experiments is 
the Fermi-LAT currently examining cosmic gamma-rays, and IceCube 0], which is already monitoring high-energy 
cosmic neutrinos. In their mandates, they both have commitments to analyze their data for the presence of dark 
matter annihilation radiation In this paper, we consider particular models from the minimal supersymmetric 

standard model (MSSM) where both neutrinos and gamma-rays are produced from annihilations; however, we also 
discuss different well-motivated models where neutrinos may be the dominant final states. These scenarios would be 
distinguishable from one another if information were received about dark matter annihilation from observations of 
both cosmic gamma-rays and neutrinos. 

Constraints on dark matter annihilation exist from observations of gamma-rays from dwarf galaxies and neutrinos 
from the Sun. It is apparent from the lack of an obvious signal from those observations that any existing indirect signal 
is faint. Successful observation of this radiation, and proper constraints from the observations, will require precise 
knowledge of other astrophysical radiation sources, and the development of rigorous predictions of the annihilation 
signal for given realistic models of particle physics and matter distribution. These predictions will be instrumental in 
understanding which kinds of models are and are not observable, and what constraints will be determined by a given 
designed experiment. In addition, this research is important for the understanding of precisely what information is 
available in a given particle model's signal, if observed by a particular experiment. This will not only aid the analyses 
of current experiments, but will also inform the design of future experiments. 

Implementing the dark matter distribution in terms of semi-analytic models allows the basic properties of the 
densest regions (the cores of lialos) to be well-represented. It is reasonable to understand the predictions produced 
by the simplest models first, learning about the physical scales most important to the signals. Additional features 
of the distribution, as seen in simulations, can be implemented in the context of these models, and their effects 
quantified — for example, different models of halo core densities, implementation of halo substructure, implementation 
of distributions of more complex halo shapes, more general halo statistics, etc. It is usefule to identify the relevant 
physical scales that determine the main properties of the annihilation signals. This is an important tool to quantify 
the robustness of the predictions to uncertainties of our knowledge of large scale structure, and tells us about the 
constraints available to be gained from the experiments. This will also guide our understanding of the annihilation 
signals calculated directly from simulation data, as in Q. 

We share preliminary results of predictions of the intensity spectrum of annihilation gamma rays and neutrinos 
produced from within and outside the galaxy for a spherical halo model distribution. In Section [Hi we describe some 
of the particle physics models that motivate our discussion. We explain the calculation of the observed spectrum of 
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gamma rays from annihilations and show results of the calculations in Sec. IIIII Of particular interest is how different 
theories can produce the same gamma-ray spectra, even though they have different annihilation modes. We show 
corresponding results for neutrinos in Sec. II VI and discuss reasons for their consideration in addition to gamma-rays, 
such as how they break the degeneracy in spectrum for different models. We provide discussion of the results in 
Sec. El 

II. PARTICLE PHYSICS MODELS 

In this paper, we discuss particle models that are supersymmetric extensions of the standard model. In these 
models, R-parity, or some similar parity property, allows only an even number of supersymmetric partner particles to 
interact on a fundamental interaction vertex. This stabilizes the lightest supersymmetric particle, which becomes the 
dark matter candidate. 



A. MSSM 

In the MSSM, the particle content is restricted to the standard model, a supersymmetric partner for each standard 
model degree of freedom, and additional Higgs fields. The supersymmetric charged higgsinos and gauginos mix to 
produce mass eigenstates called charginos, and the neutral higgsinos and gauginos produce neutralinos. In addition 
to providing a dark matter mass candidate, this model also stabilizes radiative corrections to the Higgs mass, and 
causes unification of the standard model forces at an energy called the grand unified theory (GUT) scale. 

The most general allowed parameter space for this model has over 100 free parameters. We will restrict our 
discussion to a small portion of these. In minimal supergravity (mSUGRA) or the constrained MSSM (CMSSM), 
supersymmetric masses are unified at the GUT scale with scalar masses having value mg, and TOi/2 being the mass 
of the gauginos @. Two Higgs fields each gain vacuum expectation values, the ratio of which is specified with the 
value tan/3. We will describe the properties of viable thermally produced dark matter models with vanishing soft 
supersymmetry-breaking trilinear coupling parameters Aq, and positive mass parameter /i that couples the two Higgs 
superfields in the superpotential. We will be focusing here on universal masses that are not so large as to result in a 
dark matter particle massive enough to produce significant top quarks from annihilations. 

In this three-dimensional parameter space of mo, mi/2, and tan/3, the parameter space is typically broken up into 
four main regions: the bulk region, the focus point (also known as hyperbolic) region, the co-annihilation region, and 
the funnel region. In these regions, the dark matter particle turns out to be the lightest neutralino Xi- 

In the bulk region, both mg and are relatively small. The neutralino is nearly pure bino (the gaugino which 

is the supersymmetric partner of the weak hypercharge gauge boson) , and annihilates predominantly to bottom anti- 
bottom quark pairs 66, secondarily to tau anti-tau lepton pairs t~^t~ (moreso at larger tan /?). These processes in the 
bulk region give the correct annihilation cross section to account for the relic density, if it were thermally produced. 

Generically, larger values of mp and mi/2 result in theories with larger mass dark matter that have smaller anni- 
hilation cross sections, and therefore would result in more thermally produced dark matter in the Universe than is 
observed today. However, when considered carefully, we see other parameter space does result in the correct relic 
density, due to different mechanisms Q, according to the parameter space of interest. 

The focus point region has a branch where mi/2 remains small and mo is allowed to increase. As mo does so, 
the lightest neutralino gains a larger Higgsino component, which opens up additional annihilation channels. Here, 
annihilation dominantly produces W^W~ bosons, with a small branching fraction also producing ZZ boson pairs 
for small to moderate tan/3. For large tan/3, the Higgsino component of the lightest neutralino is again small in this 
region, but Bino annihilation is enhanced by an increased coupling to the pseudoscalar Higgs A and annihilation is 
again dominated by 66 and r+r". Some parts of this parameter space have restrictions from the results of the CMS 
and ATLAS [l3| experiments at the Large Hadron Collider. 

There is a threshold where mo becomes too small and one of the supersymmetric partners of the tau (stau f) 
becomes the lightest supersymmetric particle, which is electrically charged and therefore cosmologically disallowed. 
This threshold increases with rai/2- Near this boundary, the f mass is only slightly larger than the mass, enhancing 
the co-annihilation interaction cross section between these particles. The f 's present in the early Universe co-annihilate 
with the Xi's, and reduces the neutralino density to the correct value. This parameter space is the stau- neutralino 
co-annihilation region [ll|- When ^0 > 0, there is parameter space at low mi/2 where a supersymmetric partner 
of the top quark (stop i) becomes lighter than x^- The stop-neutralino co-annihilation region [T2i| is near this 
boundary. In these parameter spaces, x? is again nearly pure bino and mostly 66 and some t^t~ are produced from 
annihilations. Because there are no f or t particles present today, they no longer contribute to annihilations and 
the effective annihilation cross section of the neutralinos is reduced from its value at freezeout. Additionally, at low 
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tan/3, annihilation is dominated by t-channel sfermion exchange, which is hehcity-suppressed [13|. The presence of a 
strong p-wave annihilation component brings the annihilation cross section up to its needed value at freezeout, but 
slow relative motions of the particles today do not allow the p-wave to contribute. In these cases, the annihilation 
cross-sections are quite small, which make the rates of annihilations low and the intensity of annihilation radiation 
much more difficult to detect. The situation becomes better at large tan /3 where annihilation via A is a stronger 
component, lifting much of the helicity suppression. 

The final parameter space, the heavy Higgs or A annihilation funnel regions, occurs where the mass of one of the 
Higgs bosons is near half the y? mass, resulting in a Breit-Wigner resonance enhancement of the annihilations at 
freezeout interaction energies [ij]. Since the resonance does not enhance the cross section today, the annihilation 
cross sections are again lower in the present epoch for dark matter models of this parameter space. 



B. Gauged U{1)b-l Model 



Another paradigm we wish to discuss, which is interesting in the context of neutrino radiation production, is the 
U{1)b-l extension of the MSSM [l^. Here, baryon number B minus lepton number L is a gauged charge with 
associated gauge boson Z' that couples to baryons and leptons, according to their B — L charges with gauge coupling 
g' . This extension requires the presence of right-handed neutrinos N"^ for anomaly cancellation, providing a natural 
framework to explain neutrino masses and oscillations. In order for this new internal symmetry to be spontaneously 
broken, we must introduce two new Higgs superficlds Tl[ and H2, standard model neutral and oppositely charged 
under B — L for anomaly cancellation. They are coupled by a new mass parameter ^' in a new term added to the 
MSSM superpotential. The physical neutrinos are light, but N'^ heavy, by the type I see-saw mechanism 
This requires a Majorana mass for the N'^, which does not obey the B — L symmetry; however, the N'^ can have a 
Yukawa coupling to another a Higgs field with lepton number —2, which we identify with i/j- This Higgs will gain a 
vacuum expectation value around 1 TeV, producing the N"^ Majorana mass and generating the appropriate neutrino 
spectrum. Thus, by defining supersymmetric partners for each of the introduced new fields and putting them in chiral 
supermultiplets, the minimal U{1)b-l extension to the MSSM has superpotential 

W = l^MSSM + 2;dN^H„L + /H^N'^N^ + /i'H'iH^ (1) 



where H,j is the Higgs superfield of the MSSM that gives mass to the up-type quarks, and L is the superfield containing 
the left-handed leptons. Note that flavor and the weak isospin SU{2)l indices have been suppressed. 

There exists parameter space in this framework where the LSP is a supersymmetric partner of N"^, the right 
sneutrino N. If the N'^ mass is less than the N mass, then annihilations could produce a large number of N'^, which 
would then decay according to the particular model considered. In any case, one would expect many direct neutrinos 
to be produced, while photons would only be produced secondarily. 

The particular model we will consider will be a parameter space where the N has a mass of 150 GeV. In this case, 
the dominant annihilation channels are the s-wave processes NN — > N'^N'^ and TV N — > N'^ N'^ via t-channel 
exchange of B ~ L neutralinos through its coupling with the gaugino Z'. The N'^, taken to have mass 135 GeV, then 
decay exclusively to 1/ and standard model Higgs h, which we took to have mass nih = 120 GeV. At this mass, the 
Higgs decays mostly to WW* bosons, and to bb quarks, each of which produce secondary photons and neutrinos. We 
will also discuss cases where N is heavier. 



III. GAMMA-RAY RESULTS 



The production of gamma-rays due to extragalactic dark matter annihilation is estimated using the halo model 
distribution of dark matter used in |18| , where each halo is specified by its mass M and observed redshift z. The 
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take typical halos to be spherical, with a Navarro-Frenk- White (NFW) density profile [2; 



Sheth-Tormen halo mass function [19| is a good approximation of the halo distributions seen in simulations. We 
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with Avir = 180 and {p){z) being the background matter density at redshift z. The distribution of scale radii 
is described by the distribution of halo concentrations c = , which is approximately described by the physically 
motivated model of [2l| over the mass scales probed by simulations. We set the minimum halo mass scale at lQ~^Mp 
Cosmological parameters used were from the WMAP7 data 



22[, and we use the linear power spectrum of [23 



0- 
to 

describe large-scale fluctuations in the matter distribution. These inputs describe the simulations well enough to give 
us a good approximation of the annihilation signal. However, as we will establish, the signal (produced from dark 
matter distributed according to our halo model) is sensitive to low mass halo properties and halo core properties, 
which are beyond the reach of current simulations, and the knowledge of which will require a good understanding of 
interactions with baryonic matter. We consider all photons emitted since the epoch of reionization, which we estimate 
to have occurred at redshift Zmax = 10. 

Given this description of the densest regions of large scale structure, the mean extragalactic intensity of gamma- 
rays from the annihilation of dark matter particles, each of mass m, at constant s-wave relative-velocity-weighted 
annihilation cross section av with annihilation spectrum is 



I^,^G{E^) = av J ^^W{{l + z)E^,z) {p^){z). 



(2) 



The intensity window function is 



WiE^,z) 
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(3) 



where t{E^, z) is the cosmic opacity to gamma-rays [2J|) and the mean square matter density is determined from 



{p^){z) = / dM^(M,z) / <l^vpl{r\M,z). 



(4) 



Note that it is possible that av depends on the relative velocity of the annihilating particles. The most prominent 
example of this, which appears in particle models, is the presence of p-wave annihilation [25| . However, this anni- 
hilation component must be very large to affect the observed annihilation signal, and thermal freezeout models for 
which this is true have very small cross sections and are difficult to see [l^. Since, in this paper, we are focusing 
on models where dark matter is thermally produced, it is reasonable to restrict our initial investigations to pure 
s-wave annihilation. However, if a special case scenario, such as resonant annihilation (26l. [27j were present, then the 
calculation needs to be modified, as described in 
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If we take our own Milky Way Galaxy dark matter halo to be a typical halo of our large scale structure model at 
mass Mq = 2 X 10^"^ Mq^ then it has scale radius rs,G = 38.0 kpc, virial radius i?vir,G = 412 kpc, and concentration 
CG = 10.8. The important parameter here for our calculation is the scale radius, since the contribution to the 
annihilation signal due to dark matter outside this radius is very small; therefore, the virial radius definition (and 
hence the value of concentration) does not significantly affect the prediction of the galactic annihilation signal. We 
estimate the solar system's position in the halo as being Rq = 8.0 kpc from the galactic center. 

With this description, the intensity of gamma-rays, due to dark matter annihilation in the galactic halo in the 
direction of angle if) from the galactic center (assumed coincident with the halo center), is typically written as 



(5) 



where the J-factor is the line of sight integration of the square dark matter density from the solar system out the 
halo m 



J(V') = / dr 
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(6) 



with 



rmax(V') = i?0 cos V' + 



i?2 



R% sin" 
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(7) 



In Fig. [U we see the contributions of the galactic and extragalactic components of annihilation to the gamma-ray 
intensity for different lines of sight in the halo. The particle physics model used in this example is the focus point region 
in mSUGRA parameter space where the dark matter particle is the lightest neutralino having mass m = 150 GeV 
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FIG. 1: The gamma ray signal from annihilating dark matter in the directions ij] — 30°, 90°, and 150° from the galactic 
center, respectively. The dark matter shown here is a 150 GeV neutralino in the focus point region of mSUGRA with 
tan/3 = 10. The dotted (blue online) line is the extragalactic component, the dot-dashed (red online) line is the 
galactic component. The solid line is the net signal. 



at tan/3 = 10, Aq = 0, and sign > 0. Here, annihilation is predominantly to W~^W~ boson pairs. The particle 
annihilation cross section and annihilation spectrum were calculated using the computer program DarkSUSY 5.0.5 
[29} . The cross section for this model is av = 1.9 x 10^^^ cm'^/s. 

In our dark matter density distribution models, the galactic component is dominant at the peak of the signal when 
looking toward the galactic center, but the contributions of the components are comparable when looking out of the 
galactic plane or away from the galactic center. It is conceivable that with slightly different choices of distribution 
parameters, the relative importance of each may be altered considerably. The relative strength of the galactic to 
extragalactic intensity at a given photon energy is 



J^,g(^^7' V') 



= / dz 



H{z){l + zfJ{^) 



-t{{1+z)E^,z) 



The important parameters then appear in the first factor of the integrand. Fig. [2] plots the extragalactic and galactic 
contributions to this factor in units of p^/Hq, where Pc is the cosmological critical density to collapse and Hq is the 
Hubble constant. The extragalactic part is relatively fiat in scale, with an area under the curve of around ASOOOp"^/ Hq. 
The convolution with the annihilation spectrum and opacity could modify the importance of this factor, depending 
on the details of those functions. One may wonder what mass scale of halos most contributes to the mean square 
density (p^). In Fig. [3l we can see that the mass integrand goes very nearly like all the way to the maximum 

mass scale, suggesting that all mass scales contribute nearly equally to the intensity. If the mass dependence of the 
Sheth-Tormen mass function correctly describes the halo distribution down to low scales, and those low-mass halos 
have density profiles well-described by NFW, with concentrations described by the model specified, then all mass 
scales are important contributors to annihilations. 

However, let us consider for the moment the effect of the neglected substructure. At Milky Way size halos, it 
is expected that substructure will increase the annihilation rate by a factor on the order of 100, depending on the 
minimum halo mass scale (soj . By definition, the smallest halos will not have any subhalos, and larger halos will have 
more and more substructure. Thus, one would expect the largest halos to contribute the most to intensity purely on 
the basis of their substructure. 

For the galactic contribution, if our galactic halo is well-described by an NFW (or similar) profile, then the value 
of scale radius Vg has a significant effect on how concentrated the dark matter is to the galactic core. Based on 
observations of stellar velocities, it is generally estimated that our galactic halo has a somewhat smaller scale radius 
than the typical radius we used (3l| . This would result in an increase in the predicted galactic intensity. 

The scaling of the density at the core is also important. On the right plot of Fig. [51 we see how the intensity formally 
diverges as the line of sight approaches the galactic center for the NFW profile. Observing a signal from toward the 
galactic center would help to better understand how the density is distributed there in our halo, and would allow us 
to test various ideas about the effects that the central black hole and baryonic cooling have on the profile. 

It is expected that the substructure observed in the simulations would increase the galactic signal by a factor of 
a few, not as significantly as the extragalactic intensity [30| . Therefore, it is not unreasonable to suppose that the 
extragalactic annihilation could dominate over the galactic signal for most lines-of-sight that are not too close to the 
galactic center. 

In summary, our estimation of the most crucial elements in these calculations, which have the greatest effects on 
the result, is: 
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FIG. 2; Left: The magnitude of extragalactic intensity is approximately proportional to the area under this curve, 
around ASOOOp^/ Hq. Right: The corresponding contribution to the galactic intensity. 




FIG. 3: The mass integrand of the mean square density at redshift z = 0, 0.5, and 1 from bottom to top, respectively. 



• the halo scale radius, the galactic value of which has an important effect on the galactic signal component, and 
the halo distribution of which affects the extragalactic signal; and 

• the inclusion of subhalos, not yet taken into account, will also increase the predicted signal, and will depend on 
the scale of minimum halo mass. 

Thus, one can conclude from this discussion that the galactic and extragalactic annihilation signals in Fig. [T] are of 
comparable intensity, due to our value of rs,G, and the lack of substructure effects. 

Unless the annihilation at the galactic core is very bright, it will be difficult to observe those dark matter annihilation 
gamma-rays originating from there because there are so many other bright sources of astrophysical gamma-rays in 
that region that have theoretical uncertainties associated with them. A less contaminated signal, for example, would 
be the consideration of the mean annihilation signal away from the core. The galactic and extragalactic components 
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FIG. 4: Left: The mean intensity for the focus point model, averaged over all directions ip > 18° away from the 
galactic center. The plot format is the same as for Fig. [T] Right: The same calculation for a 150 GeV dark matter 
particle that annihilates purely to W~^W~ , bb, or r+r^. 
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FIG. 5: The neutrino signal from annihilating dark matter in the directions tp = 30°, 90°, and 150° from the galactic 
center, for the same particle model and plot format as in Fig. [1] 
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FIG. 6: The mean neutrino intensity for the focus point model. Left: The all-sky intensity, 0° < V' < 180°. Middle: 
The anticore intensity, ip > 18°. Right: The core intensity, ip < 5°. The plot format is the same as in Fig. [TJ 



for this are shown in Fig. |4] for the same focus point model. For comparison, we also show total signals for dark 
matter, of the same mass, that annihilates to W^W^ , bb, or r+r^, at the same annihilation cross section as the 
focus point model. The annihilation spectra for these models were simulated with the event generator Pythia 6.135 
[3^ 1. The sources of photons in these models are from decaying pions or radiating charged fermions. The W and b 
spectra are more dispersed to lower energies because they are more likely to decay to hadronic showers where each 
photon-emitting product is at lower energy. At 150 GeV dark matter annihilation, the photons from annihilation to 
W^W~ are indistinguishable from annihilation to bb. These pure branching ratio intensities can be used to construct 
the intensity profile for any theory that annihilates to these states, with known branching ratios. For larger dark 
matter masses, the W and b signals become more distinguishable from one another. 



IV. NEUTRINO RESULTS 

Because the models discussed in the previous section also contribute a neutrino annihilation spectrum it is 

interesting to consider this component of the signal as well. Because the neutrino is electrically neutral and weakly 
interacting, it also propagates relatively freely through the cosmos, and the annihilation signal will have both galactic 
and extragalactic contributions. This calculation is completely analogous to that for the gamma-ray signal. We 
neglect any cosmic opacity for the neutrinos in the sample calculations that follow. 

In Fig. [5l we show the galactic, extragalactic, and net intensity of cosmic neutrinos from annihilations of the same 
150 GeV focus point neutralino dark matter considered in the previous section. In the galactic signal, we clearly see 
the peaks from primarily and secondarily produced neutrinos from the W decays. However, those features are washed 
out in the redshift-modulated extragalactic signal. We note how both components contribute significantly to the total 
signal in all of the shown lines of sight. Again, reasonable adjustments of dark matter distribution parameters and 
consideration of halo substructures could significantly alter this balance in either direction. 

Although the neutrino signal still suffers from uncertainties in the galactic core density profile, it does not suffer 
from the same astrophysical contamination as do gamma-rays. Therefore, there is no reason to exclude the galactic 
center in these experiments. In fact, if a neutrino detector with high angular resolution can be developed, it is a good 
strategy to focus on the galactic center. 

Fig. [6] shows the neutrino signal for the focus point model averaged over the whole sky, directions away from the 
core, and directions focused on the core, respectively. We see how the galactic signal is seen to dominate the signal 
at the galactic core if we assume the NFW profile holds to the center, and we neglect extragalactic substructures. 
The same dominance of the galactic core occurs with annihilation gamma-rays, but it is very difficult to see those 
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FIG. 7: All-sky neutrino plus antineutrino detection rates for 150 GeV dark matter annihilation. The thick lines are 
for electron or muon flavor, and the thin lines show the tau flavor rate. At 70 Gev, the top two lines (green online) 
are for annihilation to r+T~ leptons, the middle two lines (blue online) show annihilation to W^W~ bosons, and the 
bottom two lines (red online) are for annihilation to bb quarks. 
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FIG. 8: All-sky neutrino plus antineutrino event rates for 150 GeV sneutrino dark matter that annihilates to two 135 
GeV right-handed neutrinos (each flavor equally represented), each of which decays to a light neutrino and 120 GeV 
standard model higgs particle. 



photons from the noisy center of the galaxy. Further work, with more realistic distributions, should better elucidate 
the situation at the galactic core, and provide an understanding of the information about the dark matter distribution 
uncertainties that may be available in an observed neutrino signal. 

It is common in the literature to express neutrino signals as binned detector event rates per detector mass. If {Iu)fi 
is the mean annihilation intensity in a solid angle of observation, the event rate for a neutrino i^f of flavor f — e, 
/i, or T in an energy bin Ei < Ei, < Ei+i is 

Rvf,i — / <iE„ a^^N{E^) {I^)^{Ei,) 

f^m J Ei 

where Na is Avogadro's number, is the molar mass of the detector material, and a^^N is the neutrino-nucleon 
charged current scattering cross section [33|. Note that NA/rim is simply the nucleon number per detector mass. To 
ease conversion of the results for different detector materials, we show the results for rim = 1 g/mol. 

Fig. [7] shows the neutrino event rates for annihilation into W bosons, b quarks, or r leptons. The logarithmic GeV 
energy bin size used is 

A = log.o = 0.04. 

At these neutrino energies, the electron and muon neutrinos have indistinguishable nucleon scattering cross sections, 
which are larger than that for the tau neutrinos. Hence, the tau neutrino event rates are a little smaller. 

Since r leptons always decay to a primary neutrino, while W bosons only decay directly to leptons some of the 
time, the v production from r's is more intense. The b quarks do not produce primary neutrinos, and only have a 
lower energy neutrino spectrum from secondary chains. Thus, the flux of neutrinos from annihilations breaks the 
degeneracy between annihilation into W~^W~ and bb that occured in the gamma-ray signal. 

Another class of models that results in interesting phenomenology for dark-matter-annihilation neutrinos is the 
U{1)b-l extension of the MSSM, described in Sec. Ill Bl in the parameter spaces where the sneutrino N is the dark 
matter particle. The neutrino detector rates for one example model are shown in Fig. [51 Here, the 150 GeV N 
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FIG. 9: Neutrino plus antineutrino event rates for 150 GeV dark matter annihilating to 2 prompt neutrinos v with 
cross section av — 1.1 x 10~^^ cm'^/s (solid lines), shown with the mean atmospheric neutrino plus antineutrino rates 
at the Kamioka site during low solar activity (dotted lines). For the atmospheric neutrinos, the upper line is the 
muon flavor, and the lower line is the electron flavor. For the annihilation neutrinos, the upper line shows the rate for 
electron flavor, as well as the rate for muon flavor. The lower line shows the rates for Vt + Vr. Upper left: the mean 
neutrino rates from the whole sky. Upper right: rates when excluding the galactic core, ^ > 18°. Lower left: rates 
when focused on the galactic core, tp < 5°. Lower right: rates when focused on the inner galactic core, tp < 1°. 



annihilates into 120 GeV right-handed neutrinos N'^ that decay to a standard model neutrino and higgs boson, the 
latter of which decays mostly to &'s and r's. The N annihilation in this model does have a slight p-wave component, 
and the s-wave cross section is av = 1.1 x 10~^^ cm^/s, giving the correct thermal dark matter relic density. 

The secondary neutrinos produced from the higgs decay result in a broad, soft spectrum, whereas the neutrinos 
produced directly from N'^ decays produce a narrower peak at lower energies on the order of the mass difference 
between the iV^ and the higgs. Due to the higgs decays, there is also a gamma-ray component to the signal. 

In the case where the N dark matter is heavier (larger than twice the higgs mass), and the N'^ mass still slightly 
smaller than it, then the physical neutrino peak occurs closer to the dark matter mass energy. This will produce a 
hard spectrum with narrow peak from the primary neutrinos, and broad low-energy tail produced by the higgs decays. 

Another intriguing scenario occurs when the dark matter annihilates solely to two light neutrinos v. In the context 
of the B — L model previously described, this corresponds to the limit where the Higgs mass is small, negligible 
compared to the N mass, and the mass difference between N and N'^ is also very small. Then the spectrum of the 
produced light neutrinos is at the energy of the N ^ and the width of the spectrum is very small. This simple scenario 
results in a prominent neutrino line feature with no corresponding gamma-ray observations. At this energy scale of 
neutrino energies, the dominant astrophysical source is atmospheric neutrinos. The solid lines in Fig. [9] show the 
detector rates for annihilation of 150 GeV dark matter particles into prompt neutrinos, to each flavor equally, with a 
cross section oiav — 1.1 x 10^^^ cw? /s with our modeled dark matter distribution. The upper line shows the electron 
flavor rates and muon flavor rates. The lower line is the tau flavor detection rate. Shown is the mean all-sky signal 
(0° < ip < 180°), an anticore signal {ip > 18°), a core signal (V' < 5°), and an inner core signal (ip < 1°). The width 
of the spectral line feature is due to the velocity distribution of dark matter in the galactic halo, which is negligible 
compared to the energy resolution of viable detectors. Therefore, it is completely contained in the energy bin at the 
dark matter mass. The diffuse component is due to the redshifted extragalactic neutrinos. The dotted lines in the 
figure are the predicted mean atmospheric neutrino rates, as would be seen at the Kamioka site during minimum solar 
activity The upper line shows the Vfj, + Vfj, rates, and the lower line shows the Ve + Ve rates. 

By comparing them with the previous neutrino rate plots, we see that the typical diffuse signals are well below 
the current measured atmospheric neutrino rates. Again, the situation likely improves with the consideration of halo 
substructure, and the background can also be reduced with respect to the signal by focusing on a nearby dark-matter- 
dense region of space, as we shall discuss for the prompt neutrino production example. 
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FIG. 10; The approximate logarithmic bin size required for the spectral line detector rate bin to reach the atmospheric 
neutrino detection rate. The upper line is for + T'e and the lower line is for +1^^. Left: ijj < 180°. Center: 
i! < 5°. Right: V < 1°- 



The prominence of the peak at the galactic core shows how a neutrino detector with high angular resolution may 
extract a spectral line feature by focusing on a dense region of space. Although the signal to background ratio improves 
with small solid angles of observation, the detection rates become forbiddingly small. With better energy resolution, 
an experiment can also gain a stronger signal in the spectral line scenario. Thinner energy bins have a higher spectral 
line height. The energy bin width at E,^ = m required for the bin height to be at the corresponding atmospheric 
neutrino rate when observing in solid angle is approximately 

Ah{m,il) ~ 



where 0g(to, ^) is the flux of galactic annihilation neutrinos of energy at the dark matter particle mass m originating 
within the solid angle il, /atm('7i) is the mean intensity of atmospheric neutrinos of energy to, and /eg('7i) is the 
mean extragalactic annihilation neutrino intensity. The corresponding required logarithmic bin width is (assuming 
AE < to) 

AE 



•AnlO' 



This approximate logarithmic energy bin size is shown in Fig. [10] for ranges of the dark matter mass, and for different 
solid angles centered on the galactic center. For comparison. Fig. |9] used A = 0.04. We see the energy scales where 
the spectral line is most hindered by the atmospheric neutrinos. At high dark matter mass, the electron neutrinos 
are much easier to see, since the electron atmospheric neutrinos are much less abundant than the muon atmospheric 
neutrinos. 

The neutrino annihilation signal is also complementary to the detection of neutrinos from dark matter annihilations 
in the Sun. While the galactic and extragalactic signals depend on the dark matter's self-annihilation cross section, 
the solar annihilation signal is primarily dependent on the elastic nucleon-scattering cross section, provided that the 
annihilation rate and capture rate of dark matter particle by the Sun are in equilibrium. This fact means that it is 
still possible to probe models for which the elastic scattering cross section is too low for annihilations to be observed 
from the Sun. For example, in the context of B — L models with N dark matter, this occurs when there is a small 
mass splitting between the real and imaginary parts of N, and the dark matter becomes the lighter of the two |17| . 



V. DISCUSSION 



The possibility of indirect detection of dark matter annihilation through astrophysical observation of its products 
is an idea that improves our understanding of the particle nature of dark matter. It can already provide constraints 
on the s-wave part of its annihilation cross-section for a range of particle masses that annihilate according to some 
specified annihilation spectrum and assumptions about the dark matter distribution. These dark matter constraints 
are not available from any other kind of experiment. 

Observation of annihilation products would provide valuable information about the dark matter self-interactions 
and distribution. The precise nature of the information that would be available to be extracted is dependent on the 
details of the dark matter properties, and on the nature of the data collected. 

In this paper, we explored how the simultaneous observation of annihilation gamma-rays and neutrinos allows for 
more constrained conclusions of the dark matter properties and distribution. We gave sample calculations of the 
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galactic and extragalctic contributions to signals from different particle physics models using a simple smooth halo 
model of dark matter distribution. It is interesting to consider gamma-rays and neutrinos because the direction 
of their source is constrained, and their signal contains contributions from both galactic and extragalactic sources. 
Observing these signals would improve the ability to reconstruct the particles produced in the annihilations, and the 
dark matter mass, by breaking degeneracies that exist in one signal. For example, dark matter of mass 150 GeV 
that annihilates directly into W^W^ bosons or bb quarks were seen to produce very similar gamma-ray spectra, but 
those different cases could be distinguishable in observations of cosmic neutrinos. For larger dark matter mass, the 
gamma-ray signals from annihilations into W^W~ will be more distinguishable from annihilations into bb. 

It is also possible for theories to be dominant in one kind of signal, and therefore undetectable in other channels. 
We discussed scenarios where annihilation could be neutrino-dominant. One intriguing case where this can happen is 
if dark matter is made up of supersymmetric partners of neutrinos, as we showed in the context of a low scale B — L 
gauge theory, which can account for narrow spectral peaks in the cosmic neutrinos. 

Once information about the dark matter particle mass is obtained, this theoretically constrains the Jeans mass 
scale responsible for the dark matter halo minimum mass. According to the halo model, the extragalactic signal 
is produced nearly democratically by all present scales of halo masses, when halo substructure is neglected. Since 
substructure has the most significant effect on large halos, one would predict the extragalactic signal to be dominated 
by the substructure of the most massive halos, which in turn is dependent on the scale of minimum halo mass. 

From this analysis, we find that the magnitude of the annihilation products would provide information about the 
annihilation cross section, the concentration of subhalos within the most massive halos, and the density profile of 
their cores. The degeneracies that these quantities have on the intensity could be broken by information in a signal 
extracted verifiably from galactic sources, if enough independent observations can sufficiently constrain the dark 
matter distribution at those sources, perhaps via the angular distribution of the galactic signal to extract information 
about a universal core density profile. Neutrinos may be much more suited than gamma-rays for observing a signal 
from the dense galactic core, since there are fewer astrophysical sources of neutrinos at the relevant energies originating 
from that region. 

It is evident that interpretation of observations will require precise, realistic theoretical predictions. It is possible 
to efficiently explore the effects of various aspects of large scale structure and particle physics with the development 
of semi-analytic descriptions of the structure. Through these methods, the important scales contributing to the 
predictions can be identified, and the robustness of the calculations against the uncertainties determined. It is thus 
that we may determine the constraints possible with modern and future experiments searching for indirect signals of 
dark matter annihilation. 



The authors thank Carsten Rott for helpful comments. The work of R.A. is supported by the University of New 
Mexico Office of Research. The work of S.C. and B.D. is supported in part by DOE Grant DE-FG02-95ER40917. 



While the density cusp at the center of dark matter halos in the NFW model causes the observed dark matter 
annihilation intensity to be infinite in the direction toward the center of the halo, mean intensities over solid angles 
including the center are finite. Baryon cooling and, in larger halos, the presence of a supermassive black hole are some 
of the important effects that ultimately generate a more realistic core profile. To keep the dark matter distribution in 
this paper relatively simple, we do not attempt to model these effects, but assume the NFW profile throughout the 
halo. 

In this appendix, we share our method for accurate calculation of annihilation intensity averages I(tPm) from 
observations over solid angles centered on the galactic center, with angular radius ipM- Refering to Eqs. ©-([T]), we 
wish to evaluate 
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Appendix A: Algorithm for mean annihilation intensities across the core of NFW profiles 
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FIG. 11: Galactic coordinates used for calculating the mean intensity due to dark matter annihilation in the smooth 
component of the galactic halo. 



Let X be the distance from the solar system, in units of the galactic halo scale radius Ts.Gi along a line of sight at 
angle tp from the galactic center, and Xq be distance of the solar system from the galactic center, also in units of r^.c- 
Then 

1-COS^M -,, , f^'^'', . , r--^''^ dx 

— 2 J[Vm) = d^/'smV' ' 

Ps,g'^s,g Jo 

where 



{x^ - 2x,x cos ^ + 4) (l + y-^-2.e.cosV. + x|) 



a;max(V') = xq cos-0 + \Jc%- {xq sin-fA)^ 



expresses the halo boundary and, as before, the halo concentration is cq = Rvir,G/ 1^3,0 ■ The integrand of J in these 
coordinates is irregular in the neighborhood oi tjj = and x = Xq, precisely where the modeled density diverges at 
the halo center. 

The accurate evaluation of this expression is more easily attained when x is replaced in favor of 9, as pictured in 
Fig.HH 

. XQsimp XQsinip 



In these coordinates. 



where 



\J {x ~ Xq cos "f/;)^ + sin^ ij) 



Xq Jo -IOmW 



Ps,g'^s,g Xq Jo JbmW Vsinfl + XQsinV' 



&m{'>Ij) — sin ( — ^ sin'0 

\CG 

The inner integration is now well-defined and easy to numerically evaluate, except for when ■0 = 0, where the path 
of integration becomes degenerate, initially at the Sun having the value of tt, and instantaneously becoming when 
crossing the galactic center. Since this degenerate point is an end of the ip integration, it is sufficient for numerical 
evaluation to consider the value of the inner integration in the limit as '0 approaches 0. 
For ^p — > 0, we have 0m — ^ xqiP/cg — > 0, and the inner integral approaches 

mW \sm0 + XQijJ 



14 



For ■0 = TT, the integration path is simply of zero measure with = constant along the path. Therefore, the 
inner integration vanishes for this value of ■(/'• 
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